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Abstract 

The sum capacity on a symbol-synchronous CDMA system having processing gain A'^ and supporting K power constrained 
users is achieved by employing at most 27V — 1 sequences. Analogously, the minimum received power (energy-per-chip) on the 
symbol-synchronous CDMA system supporting K users that demand specified data rates is attained by employing at most 27V — 1 
sequences. If there are L oversized users in the system, at most 27V — L — 1 sequences are needed. 27V — 1 is the minimum number 
of sequences needed to guarantee optimal allocation for single dimensional signaling. 7V orthogonal sequences are sufficient if a 
few users (at most 7V — 1) are allowed to signal in multiple dimensions. If there are no oversized users, these split users need to 
signal only in two dimensions each. The above results are shown by proving a converse to a well-known result of Weyl on the 
interlacing eigenvalues of the sum of two Hermitian matrices, one of which is of rank 1. The converse is analogous to Mirsky's 
converse to the interlacing eigenvalues theorem for bordering matrices. 

Index Terms 

Code division multiple access (CDMA), eigenvalues, interlacing eigenvalues, inverse eigenvalue problems, sequences, tight 
frames 

I. Introduction 

Consider a symbol-synchronous code-division multiple access (CDMA) system. The fcth user is assigned an TV-sequence 
Sk G of unit energy, i.e., s\sk = 1. The processing gain is TV chips, and the number of users is K, with K > N. User fc 
modulates the vector Sk by its data symbol Xk G M and transmits XkSk over TV chips. This transmission interferes with other 
users' transmissions and is corrupted by noise. The received signal is modeled by 

K 

Y = Y, SkXk + Z, 

k=l 

where Z is a zero-mean Gaussian random vector with covariance In, the TV x TV identity matrix. 

We consider two problems already studied in the literature. Problem I: User ft has a power contraint pj. units per chip, i.e., 
E[X'^] < Npk- The goal then is to assign sequences and data rates to users so that the sum of the individual rates at which the 
users can transmit data reliably (in an asymptotic sense) is maximized. The maximum value Csum is called the sum capacity. 
Viswanath and Anantharam [1] show the following. Let ptot = SfcLiPfc- 

• Oversized users, i.e., those capable of transmitting at large powers relative to other users' power contraints, are best 
allocated non-interfering sequences; 

• others are allocated generalized Welch-bound equality (GWBE) sequences [1]; 

• Csum < (1/2) log(l +Ptot)' with equality if and only if no user is oversized; 

• no user is oversized if TVp^. < ptot for every user k. 

Problem II, a dual to Problem I, is one where user k demands reliable transmission at a minimum rate rk bits/chip. The 
goal is to assign sequences and powers to users so that despite their mutual interference and noise, each of the users can 
transmit reliably at or greater than their required rates, and the sum of the received powers (energy/chip) at the base-station is 
minimized. Guess [2] shows results analogous to those of Problem I. 

• Oversized users, i.e., those that demand large rates relative to other users' requirements, are best allocated non-interfering 
sequences; 

• others are allocated GWBE sequences; 

• the received sum power is lower bounded by exp {2rtot} — 1, where rtot = '''fe' ^^'^^ equality if and only if no user 
is oversized; 

• no user is oversized if Nrk < rtot for every user k. 

Once the optimal sequences are identified and the power or rate allocated to a user determined, the quantities have to 
be signaled to the typically geographically separated users. This "control-plane" signaling eats up some bandwidth on the 
downlink. This can be a significant fraction of the available resources when the system is dynamic: users may enter and leave 
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the system, or the channel may vary as is typical in wireless channels. The sequences and allocations may need periodic 
updates. Every such update requires the transmission to each user of an A^-vector representing the sequence and a real number 
representing the power or rate allocated. It is therefore of interest to identify schemes that result in reduced signaling overhead. 

In this paper, for both Problems I and II, we come up with a generalized WBE sequence allocation for the non-oversized 
users. The allocation employs at most 2N — L — 1 sequences, where L is the number of oversized users. To prove this, we 
only need to restrict our attention to the case when no user is oversized and show that we need at most 2N — 1 sequences. 
The 2N — L — 1 requirement in the presence of oversized users follows immediately. This upper bound on the number of 
sequences enables reduced signaling on the downlink when there are a large number of users relative to the processing gain. 

Following Guess [2], we show achievability with such an allocation via a successively canceling decoder. Indeed, if there 
are fewer sequences than users, two or more users will be assigned the same sequence and they will overlap with each other. 
Nonlinear receivers, among which the successively canceling decoder is an example, are therefore necessary. Our algorithm to 
identify the sequences and powers or rates allocated requires 0{KN) floating point operations and is numerically stable. 

The solution to the above problem draws from an interesting result in matrix theory. If A is a matrix, let A* denote its 
Hermitian conjugate. If A is Hermitian, let <j{A) denote its eigenvalues in decreasing order. The following is a well-known 
result due to Weyl (see for e.g., [3, Section 4.3]). If A is an iV x Hermitian matrix with a{A) = \ — (Ai, • • • , Xn), and 
c £ C^, then (j{A + cc*) = ^Ai, • • • , A^^ = A satisfies the interlacing property 

Ai > Ai > As > As > ■ • • > Aat > Aat. (1) 

In this paper, we demonstrate a converse to Weyl's result. We show that for any Hermitian A with a{A) = A, and any real A 
such that A and A satisfy ([T]), there exists c £ such that (t{A + cc*) = X. This converse facilitates a sequence allocation 
algorithm that, for both Problems I and II, satisfies the declared properties on the number of sequences. The above problem 
falls within the class of structured inverse eigenvalue problems. Chu [4], Chu and Golub [5] provide excellent reviews of 
similar inverse eigenvalue problems. The problem we have solved appears to be new. 

We then investigate the need for 2N — 1 sequences. While there are rate requirements (respectively, power constraints) for 
which fewer sequences suffice, in general, 2N — 1 is the minimum number required to meet any set of power (respectively, 
rate) constraints. 

Finally, we show that if some users can be split, N sequences suffice. This is useful because any set of N orthogonal 
sequences will then work. For example, we may employ the standard basis as in a TDMA system, or a set of Walsh sequences 
when is a power of 2 as in a CDMA system. Moreover, this set can be fixed up front. As the system evolves, it is sufficient 
to send an index to this set, thereby making the sequence signaling on the downlink rather easy. Some users may be signaled 
more than one index. If there are no oversized users, at most — 1 users are split into exactly two each, and therefore will 
need two indices. Because their energy is now concentrated in a two-dimensional subspace instead of one, the benefits of 
spreading (such as robustness to jamming) is obtained to a lesser degree. The allocation is however optimal; it will either 
maximise sum capacity or minimise sum power. 

It is perhaps obvious that if multi-dimensional signaling is allowed, then N orthogonal sequences suffice. Indeed, the goal of 
an optimal allocation is to ensure that energy is spread equally in all dimensions leading to a GWBE allocation. A sufficiently 
fine splitting of the users into virtual users with smaller requirements, coupled with multiple dimensional signaling per user, 
will achieve this. The interesting aspects of the above result, however, are an identification of the number of users that need 
be split and the resulting dimensionality of their signaling. 

The problem of sequence detection has attracted much attention since the work of Rupf and Massey [6] who consider the 
scenario where users have identical power constraints and the goal is to maximize the sum capacity. Viswanath and Anantharam 
[1] extend the results for users with differing power constraints. Guess [2] studies the dual problem (Problem II) and obtains 
sequences that minimize sum power subject to users' rate constraints. Waldron [7] shows that GWBE sequences are tight 
frames, a generalization of bases useful in analyzing wavelet decompositions. Tropp and others [8] show that sequence design 
for Problem I is a structured inverse singular value problem [5] and provide a numerically stable algorithm whose complexity 
is 0{KN) floating point operations. Our proposed algorithm (first presented in [9]) has the same 0{KN) complexity and 
numerical stability properties. Moreover, our algorithm guarantees optimality with at most 2A^ — 1 sequences and will work 
for any ordering of users. 

Several iterative algorithms for finding optimal sequences have been proposed. As our focus is on finite-step algorithms, 
we refer the interested reader to a work of Tropp and others [10] and references therein. Another line of research has been 
that of sequence design for suboptimal linear receivers. Viswanath and others in [11] focus on linear MMSE receivers and 
identify sequences and powers to meet a per-user signal to interference ratio (SIR) requirement. Guess [12] extends the results 
to decision-feedback receivers. In our work, we place no restriction on the receivers. 

The paper is organized as follows. Section |ll] provides the preliminaries and states the problems. Section |lll] discusses the 
converse to Weyl's result. Section |IV] provides the new algorithms, shows their optimality, and verifies that at most 2N — 1 
sequences are utilized. It also discusses the complexity and numerical stability of the algorithms. Section |V] deals with the 
need for 2N — 1 sequences and Section [VT] discusses the rate-splitting approach to minimize the number of sequences. Section 
I VIII provides some concluding remarks. 
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II. Preliminaries and Problem Statements 

Suppose user k is assigned sequence Sk and is received at power pk- Let S be the N x K matrix [si S2 ■ ■ ■ sk]- Then, the 
capacity region [13] can be written as 

C{S,p)^ n {(ri,...,rK) eM^ : 

JC{1,---,K} 



fee,/ 



■ ■SfeSfc 

fee J 



(2) 



where |A| denotes the determinant of the matrix A, and is user fc's data rate in bits/chip. 
The rate vector r = (ri, • ■ • , tk) is a verfex of this capacity region if 



^fe = ^log|^fe|-^log|^fe-i|, l<k<K, (3) 



where Aq = /^r, and A^ = + Npf^ ■ Sfcs|.) for k — 1, ■ ■ ■ , K. The matrix Aj. is real, symmetric, and positive-definite. 

The vertex (cf. ^) satisfies > and r G C{S,p). These are deduced as follows. Let ^{Ak) ~ (^X'i \---,x'}^^y Clearly 

A^°^ = • ■ • = A^^ = 1. Weyl's result O indicates 

^(fe) > ^(fe-i) > ^(fe) > ^(/c-i) > . . . > A« > Aiv'-'^ > 1, 

which implies that \Ak\ > and hence > 0. Moreover, r G C{S,p) because this rate point can be achieved via 

successive decoding. (Alternatively, C{S,p) is a polymatroid [14], and therefore contains all its vertices in ^^). 
With the above ideas fixed, let us now re-state Problems I and II. 

• Problem I : Given a per user power contraint of p = {pi, - ■ ■ ,pk) where no user is oversized (i.e., Npk < ptot for every 
user fc), find S and r that satisfy r G C{S,p) and rtot = (1/2) log(l +ptot), the maximum sum-rate among all sequence 
and rate allocations. 

• Problem II : Given a per user rate requirement of r = (ri, • • • , r^f ) bits/chip where no user is oversized [i.e., Nr^ < rtot 
for every user fc), find 5 and p that satisfy r G C{S,p) and ptot = exp {2rtot} — 1, the minimum value among all power 
and sequence allocations. 

Section |IV] identifies the matrices Ak and the sequence matrix S that solve Problems I and II. Vertices will play an important 
role in the solution. 

III. Converse TO Weyl's RESULT 

Proposition 1: {Converse to Weyl's result) Let A be an x complex Hermitian matrix with (j{A) = (Ai, • • • , Xm)- Let 
(Ai, • • ■ , Xn) be N real numbers satisfying the interlacing inequality Then, there exists a c G such that (j{A + cc*) = 
(Ai, • • • , Xm)- If, in addition, A is a real, there exists c G with the same properties. □ 

This result is a pleasing dual to an existence result shown by Mirsky [15] for bordered matrices. The problem falls under 
the category of additive inverse eigenvalue problems [4], [5]. Proposition [T] however appears to be new. A proof can be found 
in Appendix I] 

Proposition [T] is crucial to this paper It guarantees the existence of a program that takes in the matrix A, a set of prescribed 
eigenvalues A that interlace the eigenvalues of A, and outputs c — c ^A, A j such that A + cc* is the desired perturbation of 
A. Let us now focus on the nature of c in a special circumstance that will be of relevance to the sequence design problem. 

Proposition 2: Let A be an iVxA^ complex Hermitian matrix with (j{A) = A = (Ai, • • • , A^v)- Let A be unitarily diagonaUzed 
by the unitary matrix U, i.e., A = UAU*, where A — diag{A}. Fix / satisfying 1 < I < N. Let A be real numbers such 
that Xj — Xj for all j except when j — I. Furthermore, let A and A satisfy the interlacing property ([TJ. Then 

1) with c = (^VXi — Xi^ ui, where ui is the Ith column of U, we have a{A + cc*) = A; 

2) A and {A + cc*) are simultaneously diagonaUzed by U. 

□ 

Proof: This is easily verified by looking at the assignment in the proof of Proposition [T] Instead of doing this, we give 
a direct proof. Let U = [ui ■ ■ ■ ujy], where Ui G C^. Then Aui = Xiui, and u*Ui — if i ^ I. Hence for every i I, 

{A + cc*) Ui ~ Aui + c (c*u,j) = Aui + = XiUi, 

and thus Ai, z 7^ /, are — 1 eigenvalues of ^ + cc*. Furthermore, 

{A + cc*)ui = Xiui + (^Xi - A;^ ui{u*iui) = Xiui, 
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and thus A; is an eigenvalue of A + cc* . This proves the first statement. 

The eigenvectors for the two Hermitian matrices are the same and therefore the same unitary matrix (of eigenvectors) 
diagonahzes them. ■ 

IV. Sequence Assignments 
Let the rate vector r, the sequence matrix S, and the power vector p be such r G C{S,p). Then 

Ptot > exp{2rtot} - 1, (4) 

where ptot — X^f^iPfe' ^^'^ '''tot — SfLi ''fc- Indeed, this can be easily seen through the following sequence of inequalities. 
Let P — diag{p}, a K x K diagonal matrix with diagonal entries given by the elements of p. 

Ptot = ti-ace P 

= ti-ace PS*S (5) 
= ti-ace SPS' (6) 

= ^ ■ trace (In + NSPS*) - 1 
1 

= — ■ trace Ak — 1 

N 

> IAkI'/"-! (7) 

> exp{2r(ot}-l, (8) 

where (|5]l follows because the diagonal entries of 5*5* are all 1; (|6ll follows because trace AB — trace BA when the products 
AB and BA are well-defined; O follows because arithmetic mean of the eigenvalues exceeds the geometric mean of the 
eigenvalues; and ^ follows because r e C{S,p) which implies that rtot < 27v \ 

For Problem I, given a set of power constraints, (HI indicates that the maximum achievable sum capacity is upperbounded by 
i log (1 + Ptot), as observed by Viswanath and Anantharam in [1]. Analogously, for Problem II, given a set of rate constraints, 
says that the minimum sum power is lowerbounded by exp {2rtot} — 1- If these bounds are achieved, then equality in ^ 
and (O indicate that all the eigenvalues of Ak are the same, and in particular. 



Ak = (1 +Ptot) In = e-Kp{2rtot}lN, 
. well-known property of GWBE sequences ([1, Section IV]). 



A. Algorithm for Problem II 

The signal and interference matrix in the absence of any users is Aq. The matrix with all users taken into account is Ak- 
Let us now assign sequences in a sequential fashion, i.e., one user after another. 

Loosely speaking, we begin with Aq and its eigenvalues A'^'^^ = (1, • • • , 1). We then fill up dimensions with energies from 
users, one user after another, jumping to the next dimension when a dimension reaches exp {2rtot} upon addition of a user. 
Each time a user is added, say user k, we ensure that A^*^^^-' and A''^^ satisfy the interlacing property. Indeed, this is guaranteed 
precisely because there are no oversized users. Moreover, they differ in at most two eigenvalues. Upon each addition, we then 
simply appeal to Proposition [1] to get Ck for user k. The sequence is then Sk — Ck/\\ck\\ and power pk = {c\ck)/N. 

The following algorithm makes this notion precise. For simplicity, we assume that users are assigned in the increasing order 
of their indices. The order of the users is immaterial as will become apparent in the proof of correctness of the procedure. 
Our assignment will make the given set of rates the vertex of some capacity region. 

Algorithm 3: Problem II 

• Initialization: Set the following: 

a','') ^1, for < fc < i^, 1 < n < A^; 
fc ^ 1; 
n ^ 1; 
Amax exp{2rtot}; 
Aa ^ In- 

• Step 1: {All users are done) . If fc > K, stop. 

• Step 2 - Case (a) : ( Only one eigenvalue changes ) . If 

A(f-i' • exp {2Nrk} < A 

max ; (9) 
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then, set 



Go to Step 3. 

Step 2 - Case (b) : ( Only one eigenvalue changes; dimension n is to be filled up) . If 



then, set 



Now set 



^n^ <- Amax, for j = fc, 

n <— n + 1. 



(10) 

(11) 
(12) 



Go to Step 3. 

Step 2 - Case (c) : {Two eigenvalues change; dimension n is to be filled up) . If we are in this step, then 



A(f-i) •exp{2iVrfe} > A„ 



Set the following: 



for j ^k,--- ,K, 



An 



(13) 

(14) 
(15) 



(k — l) 

Observe that X)^_^_l — 1 . Now set 



n <— n + 1 . 

Step 3 : {Assign power and sequence fijr user k) . Identify the vector 

Cfe -c(Afe_i,A(^)) 

from Proposition [T]. Then set 

Sfc ^ Cfc/||cfe||, 

Pfc <- (4cfc)/^- 

A sequence and power are now allocated to user k. Finally, 

Ak ^ Ak-i+Ckcl, 
k ^ k + 1. 

Go to Step 1. 



□ 



We now make some remarks on the computations and the numerical stability of the algorithm. The matrices A^ need not 
be explicitly computed. It is sufficient to identify and store the unitary matrices Uk ~ ul ■ ■ ■ u)^ that diagonalize Ak- 

Computation of c (^Ak-i, A'''') utilizes only Uk-i, X'^'^'^^K and the new A''''\ as seen in the proof of Proposition [T] For Step 
2(c) at most two eigenvalues change in going from Ak-i to Ak- So the two matrices share N — 2 eigenvectors and exactly 
two eigenvectors change. These are computed via a rotation in the (n, n + l)st plane as follows: 



a -f3 
(3 a 



where 



a 



An ^ A„+i ) I A„ — A„ 



^ {Xn — A„_|_i^ (A„ — A„+i) 



fj = \J\ - a2 = 



An+i — A„_|_i I A„ — A. 



^ ^A„ — A„_(-i^ (A„ — A„+i) 



with A ~ and A = A''^-' for simplicity. Moreover, pk = A„+i + A„ — (A„+i + A„), and 

t-fc — (Jn U-n ^ i/n+1 "n+1 i 
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where 



An — A^ 



Vn = 



\ 



An — An 



|A„ — A„+i I 



Vn 



\ 



An+1 — A„ 



An+1 ^ A- 



|A„ — A„+i I 



These facts can be verified by direct substitution. Furthermore several assignments in (fT4l i and in the initiahzation step can be 
imphcitly assumed. The upshot is that computation of sequences as well as the new eigenvectors requires only 0{N) floating 
point operations. Steps 2(a) and 2(6) do not result in a change in eigenvectors. The sequence in these simpler cases is a scaled 
version of ul*^ ^' and requires only N floating point operations for computation. Hence the algorithm's complexity is 0{KN) 
floating point operations. 

The rotation requires identification of a and (3; these are between and 1 and hence their computation is numerically stable. 
Furthermore, it is easy to verify from the interlacing inequality that 



{( 



An ^ An 



An 



A 



n+l 



< A„, — Al- 



and therefore computations of ?/„ and ?/„+i are also numerically stable. This verifies that the algorithm is numerically stable. 
Before we give the proof of correctness, we illustrate the algorithm with an example. 



Example 4: Let N = 2. Consider four users with rate requirements ri, r2, ra, and r4, such that rk < rtot/N, for fc = 1, 
i.e., there are no oversized users. 



■A, 



Fig. [U represents the eigenvalue space of the matrices Ak, fc = 0, 1, • • ■ , 4. Our interest is in the upper sector in the positive 
quadrant which represents Ai > A2 > 1. The algorithm results in a sequence of (^\^^\x^^^, the eigenvalues of A^- 
The given rates should be supportable. This places a condition on the product of the eigenvalues 



> 



exp 



, fc = 1,2,3,4, 



i.e., after user k is added, the eigenvalue pair should lie beyond the kth hyperbola. Observe that N + Nptot = trace Ak ~ 



A^'^-' + A2"''' is least, and therefore ptot is least, if A^"''' = A^"''' = Amax = exp{2rtot} as is indicated by the tangent to the 
outermost hyperbola in the figure. An allocation that reaches this point is a GWBE allocation. The proof below shows that 
this is possible if none of the users are oversized. 

As the users are added sequentially, the interlacing condition says that the eigenvalues should lie in infinite half-strips shown 
in Fig. [T] The condition for A: = 1 is degenerate; it is the half line represented by A2 = 0, Ai > 1. 

The first user is added with minimum power to meet the interlacing and rate conditions. The other users are added to minimize 
the number of non-unit eigenvalues (as is the first user). A new eigenvalue is enlarged only if the previous eigenvalue has 
reached A„iax- The choices of eigenvalue pairs that follow this rule are indicated in the figure as solid circles. The steps 
executed for the four users are 2(a), 2(c), 2(a), 2(6), respectively, for the four users. 

Note that the choice is not greedy in the powers that are allocated. Such a "power-hungry" allocation is not guaranteed to 
reach the minimum power point, whereas the above algorithm guarantees such an allocation. □ 

Proof: We now give the proof of correctness of Algorithm [3] When a new user, say user k, is added, two conditions on 
the set of eigenvalues of the signal and interference matrix should be satisfied. One, the added user's rate should be supportable 
and hence. 



AK) _ .{K) _ 



1 

2N 



1 



N 



— logTTA, 

(=1 



(fe) 



this imposes a lower bound on the product of the eigenvalues. It is easy to verify that the algorithm satisfies the above inequality 
with equality (see JTSl)). Two, the interlacing eigenvalues condition should hold because the signal and interference matrix Ak 
is an as yet unknown rank-1 perturbation of the matrix A^^i. 



Interlacing eigenvalues : We first ensure that an execution of Step 2 results in a set of eigenvalues A^*^^ — ^A^'^^ , 
that interlace the eigenvalues A^'^^^-' = ^A^'"~^^ , • ' ' : A^ For a particular k, exactly one of the cases (a), (6), or (c) holds 
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Consider case (a) or (b). Then exactly one eigenvalue changes between A'^'^ and A^*^-*. Indeed, either (|9j or (fTTT i holds, and 

xC^-i) - i \ ... X xC^-i) 1 1 

V n-1 N~ 

A('=) - I A • • • A A^'^) 1 • • ■ 1 

y . ' ^ — ^ 

To check that the interlacing inequality holds, it is sufficient to check that 

A,„ax > A(f) > Al'^-i) > 1. 

But this is easily verified. Inequality (a) holds because of (|9]l and (fTOl i in case (a) and because of (fTTT l and ( fT2] l in case (6). 
Next, inequality (/3) holds because > in both ( fTOl i and ( fTSl i. Furthermore, inequality (7) is preserved at all steps. 

To check that the interlacing inequality holds under Step 2 case (c), note that if this step is executed, (fTST l holds, and the 
eigenvalues are given by 

aC^-i) - I A ... A A^'^-^) 1 1 ... 1 1 

AC^') - I A ■■■ \ \ X'^'"'' 1 . . 1 

^ — I ^max: : ^max; ^maxi ^7 : ^ 

\ n-1 N-n-1, 

To see that the interlacing inequality holds, it is sufficient to check that 




A. 



n+l 



This is can be verified as follows, (a) holds because of (p^ . At all times, our choice of eigenvalues is such that 1 < aI'^ < 
Amax and therefore (/3) holds. To check (7) observe that ( fTSI l gives 

,W _ Aiy^-Ar^)-exp{2jVr,} ^^(,_,^ 

^n+l — A, ^ 

'^max 

because A^j^i^^ — 1 and exp{2iVrfc} < Amax = exp{2rfot} for there are no oversized users. Lastly (6) follows because of 
^ and (O. 

Termination : After an execution of Step 3, k is the index of the next user to be allocated, and it is easy to check via 
induction and ([TOll, or dTT} and ([T2]i, or ([Tl]i and ([15), that 



A-il =exp<^ 2iV> r, V. (16) 




This ensures that every one of the K users can be added without the eigenvalues exceeding Amax at any stage and in any of 
the dimensions because 

\Ak-i\ = exp |2A^5Z I - 

for every k. In particular, if at the end of some Step 3 execution we have n — N for some k < K, then users k, k + 1, ■ ■ ■ , K 
can be added without the need for an additional dimension, i.e.. Step 2(c) will no longer be executed. Moreover, because of 
( |T6] l. we have \Ak \ — exp {2Nrtot} and therefore the last Step 2 to be executed will be Step 2(c) thereby correctly filling up 
the last eigenvalue at termination. 

Validity of the sequential procedure : The previous steps indicate that the interlacing inequality holds at every step of the 
algorithm, and therefore a sequence matrix S and a power allocation p will be put out by the algorithm. The given r, by 
design, is a vertex of the capacity region C{S,p). As C{S,p) is a polymatroid, it contains all its vertices, and therefore 
r e C{S,p). (Alternatively, a successive cancellation receiver indicates that r is achievable). This implies that the allocations 
S and p support the given rate vector r, and the proof of correctness is complete. ■ 
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B. Algorithm for Problem I 

We now look at the dual problem of sequence and rate allocation for sum capacity maximization given users' power 
constraints p — (pi, • • • ,pk)- Guess [2] explores the duality between these two problems. 

Algorithm 5: Problem I 

• Initialization: Set the following: 





- 1, 


k ^ 


- 1; 


n ^ 


- 1; 




- 14 


Ao ^ 


- In 



• Step 1: (All users are done) . If fc > K, stop. 

• Step 2 - Case (a) : ( Only one eigenvalue changes ) . If 

Xi^-'^ + Npk < Xm.^, (17) 

then, set 

Ai'^ ^Al'^-i^+iVpfe. (18) 

Go to Step 3. 

• Step 2 - Case (b) : ( Only one eigenvalue changes; dimension n is to be filled up) . If 

Al'=-i)+A^Pfc = A„,ax, (19) 

then, set 

A^f'^ ^A,„ax for j = k,---,K. (20) 

Now set 

71 <— n + 1. 

Go to Step 3. 

• Step 2 - Case (c) : [Two eigenvalues change; dimension n is to be filled up) . If we are in this step, then 

Xi^-''^ + Npk > An,ax. (21) 

Set the following: 

Al^') ^ A,nax for j = k,---,K; (22) 
A^i ^ Ai';Y^+Ai'^-i)+iVpfe-A„,ax. (23) 

Observe that A,^f_|r/^ = 1 . Set 

n ^ n + I. 

• Step 3 : (Assign rate and sequence for user k) . Identify the vector 

Cfe = c(^fc_i,A('=)) 

from Proposition [T|. Then set 

Ak ^ Ak-i+Ckcl, 

Sk Cfc/||cfc||, 

'^'^ ^ ^log|Afc| - — log|Afc_i|. 

A sequence and rate are now allocated to user k. Finally, 

k^k+1. 

Go to Step 1. □ 
We now illustrate the algorithm with an example. 

Example 6: Let N — 2. Consider a numerical example with four users having power Umits p = {pi,P2,P3,P4) = (2, 2, 3, 1). 
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It is easy to check that no user is oversized. Refer to Fig. |2] the eigenvalue space of the matrices Ak for fc = 0, 1, • • ■ , 4. 
Because trace = N + N Pj^ 'he eigenvalues are constrained to lie within the confines of the four triangular regions 
bounded by 



k 



Ai > 1, A2 > 1, Xi+ X2 < N + N2_^pj 
as given in the figure. The last bounding line Ai + A2 = + Nptot is a tangent to the region 

A1A2 > {l+Ptotf 

at the point Ai = A2 = 1 + ptot- If this point can be reached, then we have a GWBE allocation and the sum rate is 

i log(l +ptot) = ^log9, 

the maximum among all allocations meeting the power constraints given by the aforementioned tangent line. 

The algorithm ensures that at each step the eigenvalues interlace and the power constraints met. As before the number of 
non-unit eigenvalues are kept at a minimum at each step. The sequence of eigenvalue pairs are 

(A«,a(^)) = (5,1), 
(A(^\Af) = (9,1), 
(A(^\Af ) = (9,7), and 
(a(^\a(^)) = (9,9) 

The interesting aspect of this numerical example is that users 1 and 2 share a sequence, say si, and users 3 and 4 share 
another sequence, say S3, and moreover, 5*53 = 0. This completely separates users 1 and 2 from users 3 and 4, attains the 
minimum ilog(l + ptot) bound, and yet gives the benefits of spreading. Any pair of orthogonal vectors suffice. We explore 
this possibility further in Section [Vl] and show how to obtain orthogonal sequences under mild rate splitting. □ 

Proof: We now give the proof of correctness of Algorithm |5] The proof parallels the proof of correctness of Algorithm 
[3j We provide it here only to highlight the steps that rely on the absence of oversized users. 
When a new user, say user fc, is added, the power constraint 

N k 

trace = ^ A^ ^ <Y.P^^ 

and the interlacing eigenvalue condition are met. The power constraint is in fact met with equality. The rate assigned to this 
user is given by We now verify these statements. 

Interlacing eigenvalues : An execution of Step 2 results in A^'"' that interlace A^'^'^^^. Indeed, for a particular k, exactly one 
of cases (a), (&), or (c) is executed. Consider case (a) or (b). Then exactly one eigenvalue changes between A'*^^^^ and A''^^ 
and 




XC^) - I \ ... \ XC^) 1 ... 1 

^ I -^max? ; ''max? ; ; ^ 

. ' 

\ n-l N—n 

To check that the interlacing inequality holds, it is sufficient to check that 

A > a''") > A^''"^) > 1 



But this is easily verified. Inequality (a) holds because of ([TtI i and ( fTSl l, or fT% and ( |20] i, as the case may be. Inequality (/3) 
holds because in (fTsT l or in (fT9] l we have pk > 0. The algorithm is such that (7) is preserved at all steps. 
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We now verify the interlacing inequality holds under Step 2(c). As a consequence of ( 1211 1. we have 

A^^-i) - I A ••■A A'*^-!' 1 1 ••• 11 



Ji-1 N-n-1. 

aC^) - I A ■■■ \ \ A^''' 1 . 1 

— I /*maxi 1 /^maxi ^maxi ^n+li 7 i 

To see that the interlacing inequality holds, it is sufficient to check that 

A.„ax Ai^-) i A(-^) AS, 1. 

This can be verified as follows. Equality (a) holds because of the assignment in ( |22] | for j — k. At all times, our choice of 
eigenvalues is such that 1 < A,^ < Amax and therefore holds. To check (7) observe that ( |23] | gives 

because a'^'^j^'^^ = 1 and 1 + Npk < Amax = 1 +Ptot for there are no oversized users. Lastly (S) follows because of ( 1211 1. the 
assignment in (l23T l. and the equality xl^_f_i^ = 1. 

Termination : After an execution of Step 3, user k is the next user to be allocated. As before, it is easy to check via induction, 
and ( fTSl ). or ( fT9] ) and ( l20b . or ( |22] | and ( l23T l. as the case may be, that 

ti-ace = N + Ny^^pj, 

i.e., the power constraints of the just added user is respected. This also ensures that every one of the K users can be added 
because 

fc-i 

trace Ak-i ^ N + nY^Pj < N + Nptot = NX^^^ 

for every k. In particular, if at the end of some Step 3 we have n = N for some k < K, then users k,k + 1, ■ ■ ■ , K can be 
added without the need for an additional dimension, i.e.. Step 2(c) will no longer be executed. 

Validity of the sequential procedure : The previous steps indicate that the interlacing inequality holds at every step of the 
algorithm, and therefore a sequence matrix S and a rate allocation r will be put out by the algorithm. The obtained r, by design, 
is a vertex of the capacity region C{S,p). C{S,p) is a polymatroid; it contains all its vertices, and therefore r G C{S,p). ■ 

V. Number of sequences 
In this section, we study the number of sequences needed to design systems that solve Problems I and II. 

A. 2N — 1 sequences suffice 

We first show that 2N — 1 sequences suffice regardless of the number of users K. If K > 2N — 1, then there are fewer 
sequences assigned than users and therefore some users share the same sequence; they will completely overlap with one other. 
But a successive cancelation receiver enables us to receive data from all such users if powers or rates are suitably assigned. 

Theorem 7: There is an optimal sequence allocation with at most 2N — 1 distinct sequences for both Problems I and II. 
Furthermore, when there are L oversized users, there is an optimal sequence allocation with at most 2N — L — 1 distinct 
sequences. □ 

Proof: Suppose that there are no oversized users. We will show that Algorithms [5] and |5] result in a design with the 
desired property. In both algorithms, a new sequence is put out if either 

• Step 2(c) is executed, i.e., a "break-out" in a new dimension n + 1 occurs, or 

• if one of Step 2(a) or Step 2(6) is executed for the first time in a new dimension n, i.e., n was updated when the previous 
user was added. 

The former of these conditions happens at most iV — 1 times because it results in a break-out in a new dimension and 
there are at most N dimensions. The latter happens at most N times because in any dimension there is only one first step. 
All subsequent steps in the same dimension do not result in a new sequence because of Proposition |2l This ensures that the 
algorithm puts out a sequence matrix with at most 2N — 1 sequences. 
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Suppose now that there are L oversized users. These users should be assigned L sequences, orthogonal to each other, and 
orthogonal to the space spanned by the sequences assigned to the remaining users. There are N — L dimensions available for 
the remaining K — L nonoversized users. An appeal to the case with no oversized users results in at most 2{N — L) — 1 for 
these K — L users. In total, therefore, we need at most L + 2{N — L) — 1 = 2N — L — 1 sequences. ■ 



B. 2N — 1 sequences are necessary for one -dimensional signaling 
We now provide a converse to the results in the previous subsection. 

Theorem 8: For Problem I, given N, there exist K > N and power constraints p — {pi, - ■ ■ ,Pk) with no user oversized, 
such that any rate and sequence allocation with 2N — 2 or fewer sequences results in a sum capacity 

not < ^log(i +Ptot) ■ 

Similarly, for Problem II, given N, there exist K > N and rate requirements r = (ri, ■ • ■ , r^) with no user oversized, such 
that any power and sequence allocation with 2N — 2 or fewer sequences results in a minimum power 

Ptot > exp{2rtot} - 1. 

□ 

Proof: Consider Problem I. Consider K — 2N — 1 users with a symmetric power constraint 

/ \ 

Ptot Ptot Ptot 



V 2N-1 1 

Let [si, S2, • • • , sk\ be any sequence assignment with 2N — 2 or fewer sequences, and (ri, • • • , rx) a rate assignment. At least 
two users share a sequence. Without loss of generality we may assume that these are users 1 and 2, i.e., si — S2- The basic 
idea is to show that users 1 and 2 put together make an oversized compound user with power constraint pi + 152 in a 2N — 2 
user system. Indeed, 

'^Ptot '2ptot ^ Ptot 

which makes the compound user oversized. The sum capacity of the system is therefore strictly smaller than i log (1 + ptot)- 
We now formalize this idea. 
A2 is given by 

A2=In+ Npisis\ + Np2S2s\ = In+ '^^f!°' sis\, 

with eigenvalues 

\ Af-l 

and therefore 

, (if) > (2) _ 2Nptot _ 2Nptot 

\ > Ai - 1 + ^ - 1 + _^ > 1 + Ptot, 

where the first inequality is a consequence of Weyl's result applied to Ak — A2 + X^fLa ^PkSks\. 
Consequently, we have the following sequence of inequalities: 

^log(l+Pt.) = ilog(^^) (24) 

> ^loglA^I (25) 

> not, (26) 

where ( l24b holds because the power constraint implies trace Ak = N + Nptot, (l25T l holds because the arithmetic mean of the 
eigenvalues is strictly larger than their geometric mean, a consequence of the fact that one of the eigenvalues A^^-* is strictly 



larger than the arithmetic mean of the eigenvalues, and ( |26] l is a necessary condition to support the rate vector r. This proves 
the first part of the theorem. 
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The proof for Problem II is very similar. Once again, we look K ^ 2N — 1 users with symmetric rate requirements 
J, _ ... J ^if^)' '^^'^ make users 1 and 2 share a sequence. Then users 1 and 2 put together form an oversized compound 

user because 

2r(ot ■^Ttot ^ '''tot 



Now observe that 



N-l 



and therefore to support the rate of the two users we need 



which leads to 

A^"^^ > A^'^ >exp{2rtot}. 

The same sequence of inequalities, (l24l i. ( |25] ). and ( |26] ), once again hold. The strict inequality in ( l25l l holds because the largest 
eigenvalue is now strictly bigger than the geometric mean of the eigenvalues. Thus ptot > Gxp {2rtot} — 1, which completes 
the proof of the second part of the theorem. ■ 

VI. Two Dimensional signaling and sufficiency of N orthogonal sequences 

Our goal in this section is to show how to achieve sum capacity (respectively minimum sum power) by using at most N 
orthogonal sequences. Thus far, each user has been confined to signal along a single dimension. In some cases, as in Example 
|6l it is possible to get an optimal assignment with N orthogonal sequences. In general we need 2N — 1 sequences to achieve 
sum capacity or minimum power. However, by letting a few users' signals span two dimensions instead of one, it is possible 
to achieve optimality with N orthogonal sequences. We start with the following definition. 

Definition 9: A vector x = {xi,X2, - ' ' ,Xfc) has a symmetric sum partition of size N, if there is a partition of the users 
{1, 2, • • • , K} into N subsets 5*1, S2, • • • , Sn, such that 

1 ^ X 
fc6S„ fe=l 

for n=l, 2, - ■ ■ ,N. The subsets 5*1, S'2, • ■ • , Sn will be referred to as the symmetric sum partition. 

Remark : If x has a symmetric sum partition of size N, no user is oversized. This is because, any user k' belongs to Sn, for 
some n, and (|27] | implies Xk' < J2k£S ^k—Xtot/N . Note that the power constraint vector p of Example |6] has a symmetric 
sum partition of size 2. 

Proposition 10: If the rate vector r has a symmetric sum partition of size N, then N orthogonal sequences are sufficient 
to attain the minimum sum power ptot = exp {2rtot} — 1- Analogously, if the power constraint vector p has a symmetric sum 
partition of size N, then N orthogonal sequences are sufficient to attain the sum capacity rtot—-^ log (1 + Ptot)- D 

Proof: We will prove the proposition for Problem II. A similar argument holds for Problem I. Let Si, 82, - ' ' j Sn be the 
symmetric sum partition of the users. An execution of the Algorithm |3] that assigns sequences and powers to all users in a 
subset Sn before assigning to users in another subset will result in an orthogonal allocation. This is because 



^2NY,rk\ 

{ fces„ J 



exp <2N } rk> = exp {2rtot} = A 



max; 



which implies that when all users in 5„ are assigned, exactly one dimension is completely filled to Amax- Users in one subset 
will therefore be assigned the same sequence, and different subsets are assigned orthogonal sequences. We will however prove 
the proposition directly. 

Let Si, 82, - ■ ■ , Sn be the symmetric sum partition of the users. Assign sequences and powers as follows: if fc G Sn, then 



Sk = e„, (28) 

'k{n 
N 



Pk - ^^[exp{27Vrfe}-l], (29) 



where 



3k (n) = exp { 2N ''j } (^0) 

j-jes,,,j<k 
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is the interference suffered by user k due to presence of other users in the same subset S'„ with a smaller index. 
It is easy to see that 

r.~loJl + ^) (31) 



27V 

is achievable via successive interference cancelation, where the highest index user in this subset is decoded first. Users in other 
subsets do not cause interference to users in this subset. Observe that the total power assigned to users in any subset 5„ is 
given by 

^ Pfc=^(exp{2rtot}-l), (32) 

where ( [32] i follows by substitution of ( |29] l and dSOl l in the left side of ( [32] ) and by observing that the resulting sum over Sn 
has only two terms that survive. 

From (|32] | we see that total power allocated to all users is 

N 

Ptot='^ ^ Pk^ exp{2rtot} - 1, 

n=l fces„ 

thus showing that N orthogonal sequences are optimal. ■ 

Having derived a sufficient condition for optimality of N orthogonal sequences, let us see how to manufacture this condition 
from a given set of power constraints or rate requirements. 

Proposition 11: Every strictly positive vector x representing power constraints or rate requirements for K non-oversized 
users can be cast into a vector x' for K' virtual users, where K<K'<K + N — \, and x' is such that it has a symmetric 
sum partition of size N . Moreover x' is obtained by splitting K' — K users into exactly two virtual users each. □ 

Proof: Consider the cumulative requirement 

Xk=y^^xi. 

i<k 

Xk is a strictly increasing function of k and Xk — Xtot- Hence there exist — 1 distinct users with indices kj, j — 
1,2, • • • ,iV - 1, such that 



fc,--i 



Xi < 



JXtot 



ki 



N 
1 



= E-'^^- (33) 



N 

i=l 



If Strict inequality holds in ( l33T l. split user fc^'s rate as 

/ jXtot 



-k, = ('^-Xk,-i) + (xk,-'^)- (34) 



If equality holds in ( l33T l leave the user as is. For users that will be split, obtain a; from x by replacing Xk^, the requirement 
for user kj, by requirements a;^. and a;^. for two virtual users, where and Xf,. are as in ( [34] l. 

It follows that x is a vector of size K , where K < K < K + N — 1, and x has a symmetric sum partition of size N. ■ 

Users whose rates or power constraints are split are assigned two orthogonal sequences and will span a two dimensional 
subspace. The design subsequently results in N orthogonal sequences. 

It is immediate that even for the case with oversized users, the oversized users will be split into at most N virtual users. 
Non-oversized users will be split into at most two virtual users. The resulting vector r has at most K + N — 1 elements and 
has a symmetric sum partition of size N. Thus N orthogonal sequences are sufficient if rate or power splitting and multi- 
dimensional signaling is allowed for some users. We make this precise in the following Proposition. The proof is identical to 
the proof of Proposition [TT] and therefore omitted. 

Proposition 12: Every strictly positive vector x representing power constraints or rate requirements for K users, regardless 
of the presence of oversized users, can be cast into a vector x of size K , where x' has a symmetric sum partition of size N, 
and K < K' < K + N -I. □ 
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VII. Concluding Remarks 

The work in this paper was primarily motivated by a desire to reduce the amount of signaling necessary to communicate 
the sequences to the geographically separated users. This signaling usually eats up precious bandwidth in the downlink. This 
is especially a problem when the channel changes, or when users enter or leave the system. These significant events trigger a 
communication of a new set of sequences and rates/powers in the downlink. 

We showed that 2N — 1 sequences are sufficient to attain maximum sum capacity or minimum sum power This results in 
some savings on the downlink when the number of users are large. The base station may broadcast the 2N — 1 sequences as 
common information and indicate to each user only the index from the set and the rate or power allocated. If there are oversized 
users, fewer sequences are needed. The allocations are optimal and the algorithms have a simple geometric interpretation. The 
algorithms take only 0{KN) steps, have the same complexity as that of [8, Algorithm 4], and the computations are numerically 
stable. We also argued that there exist rate or power constraints that necessitate 2N — 1 sequences. 

We then saw that with a small penalty in the spreading factor for a few users, N orthogonal sequences are sufficient. The 
users and the base station can then agree on a fixed orthogonal set of N sequences. The base station only needs to signal the 
powers/rates and the indices corresponding to the sequences allocated to a user. 

The proposed algorithms identify a tight frame that meets the supplied constraints. They may therefore find application in 
other areas where tight frames arise. 

The assumed multiple access channel model, however, has severe limitations. The uplink wireless channel typically suffers 
from multipath effects, fading, and asynchronism. Moreover, the users are not active all the time. Yet, if the users can all be 
synchronized via, for example, a Global Positioning System (GPS) receiver, our results give some interesting design insights. 
For frequency-flat slow fading channels where all the users with a tight delay constraint have to be served simultaneously, multi- 
dimensional signaling (more commonly referred to multi-code where the code is a spreading sequence) can allow communication 
at sum capacity or minimum power. Orthogonal sequences are sufficient and the signaling in the downlink is significantly 
reduced. A successively canceling decoder is necessary, but the complexity of this receiver is reduced to a great extent because 
optimal decoding for the K users decouples into decoding for N separate non-interfering groups. 

Fairness of the allocation has not been considered in this paper. However, with the rate splitting approach that separates 
virtual users into groups, the decoding order of virtual users within a group can be cycled to get a fairer allocation. The first 
user to be decoded in a group treats all others in the group as interference and suffers the most. Cycling ensures that this and 
other such disadvantageous positions are shared in time by all users. 



Appendix I 
Proof of Proposition[T] 

Proof: Observe that the complex Hermitian matrix A can be unitarily diagonalized as follows 

A = UKU*, 

where U = [ui U2 ■ ■ ■ uat] is the matrix of eigenvectors of A and A = diag{Ai, • • • , Xn}, a diagonal matrix with real entries. 
Assuming that the Proposition holds for real and diagonal matrices, we can find a real vector y such that a{A + yy*) — 
(^Xi,---,XNy Clearly, 

a{U{A + yy*)U*)=<y{A + yy*), 

and therefore the desired rank-1 perturbation of A with the prescribed eigenvalues is A+ (Uy) (Uy)*, i.e., c = Uy. If A is a 
real and symmetric matrix, U may be taken to be a real orthogonal matrix, and therefore c is a real vector 

It is therefore sufficient to prove the Proposition for diagonal matrices with real entries. Let A now be such a matrix. Let 
y £ M.-^ be a nonzero vector. Let yN 0- The characteristic polynomial p^{t) of A — A + yy'^ is computed as 

tiN - A 

t- ^i-yl -ym ■ ■ ■ -ywN 
(a) -2/22/1 t-X2~yl ■■■ -y2yN 

-yNyi -yNy2 ■■■ t- Xn -ylf 





-yNV 


(«-Aiv)^.T 

Vn 


t- Xn ~y% 
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A 



N-l) 





-VNV 


(t-> 
yi\ 

In-1 


I 

VN{tlN 


-1 ~AN-iy^v 


0^ 


1 



id) 



Hn-I — Aat-i 





0^ 


(i-AA.)(l-EL?/fdl-) 



2—1 \ 2 — 1 / 



(35) 



where (&) follows by setting the columns of the matrix in (a) to be hi,i = 1,---,N, then by replacing column hi by 
hi - hNUi/vN, for i = 1, • • • ,iV - 1, and then by setting v = {yi, ■ ■ • ,?/jv-i)^ e R^^^ and Ajv-i = diag(Ai, • • • , Aat.i); 
(c) follows by multiplying the matrix in (&) on either side by triangular matrices that have determinant 1; (d) follows by 
multiplying out the matrices in (c). Observe that the conclusion remains unchanged if y^ = 0, but instead some ^ for 
some other k; such a k exists because y is nonzero. 
Define the A^th degree polynomials 

fit) ^Uliit'\), 

m =n.=i(^-A.). 

Suppose now that the eigenvalues in A are all distinct. By the Euclidean algorithm, we have 

f{t)=g{t)+r{t), 

where r{t) is a polynomial of degree at most — 1. We also have /(A^) — r{Xi) for i = 1,2, ■ ■ ■ , N because .g(Aj) — 0. 
The polynomial r{t) being known at different points can be written explicitly by using the Lagrange interpolation formula 

N .X 



Thus, 



m 



9{t) git) ^ g'iXi) t - A, 



N 



(36) 



Comparing dSST l and ( |36] |, it is clear that we are done if we can set yf — —f{Xi) / g'{Xi), for i — 1, • ■ • , iV. To do this, however, 
we must show that fiXi)/g'{Xi) < 0. Interlacing property ensures that this is indeed the case. Indeed, we can easily see that 

N 

fix.) = (-i)^-'+^ni^*-^^-i 

N 

and therefore f{Xi) and g'iXi) are always of opposite signs. Thus we may take 



n;LiiA.-A,i 

\ Uf=i,j^i |A« - Ajl 



(37) 
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Next suppose that A has muhipUcities. Observe that A and A interlace. Consequently, if the spectrum A has multiplicity k 
for a particular value 9, i.e., Xi+j ~ 9 for some / and for j — 1, - ■ ■ ,k, then the spectrum A has multiplicity fc — 1, or k, or 
fc + 1 for that 6. This is because — 1, • ■ • , fc — 1, are pegged at 6 due to the interlacing property. So the multiplicity of 

6 is at least fc — 1. It may further happen that Xi — 9 or Xi+k = or both, leading to the other possibilities. 

We can therefore set 

where the polynomials / and g have the same degree L, with 1 < L < N, and have no common factors. Moreover, g has no 
multiplicities and the zeros of / and g interlace. The same arguments that lead to ( l37l ) hold with yi^ set as in dJTl ) for those 
indices that survive and yj = when j ^ li,i = 1, ■ ■ ■ , L. The details are easy to fill and therefore omitted. ■ 
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Fig. 1. Rate constrained sequence allocation algorithm 




Fig. 2. Power constrained sequence allocation algorithm 



